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SUMMARY 


An analytical investigation is made of a precession- type instability 
which can occur in a flexibly supported aircraft-engine — propeller com- 
bination. By means of an idealized mathematical model which is comprised 
of a rigid power-plant system flexibly mounted in pitch and yaw to a fixed 
backup structure, the conditions required for neutral stability are 
determined. The paper also examines the sensitivity of the stability 
boundaries to changes in such parameters as stiffness, damping, and asym- 
metries in the engine mount, propeller speed, airspeed, Mach number, pro- 
peller thrust, and location of pitch and yaw axes. Stability is found 
to depend strongly on the damping and stiffness in the system. 

With the use of nondimens ional charts theoretical stability bound- 
aries are compared with experimental results obtained in wind-tunnel 
tests of an aeroelastic airplane model. In general, the theoretical 
results, which do not account for wing response, show the same trends 
as observed experimentally; however, for a given set of conditions cal- 
culated airspeeds for neutral stability are consistently lower than the 
measured values. Evidently, this result is due to the fact that wing 
response tends to add damping to the system. 


INTRODUCTION 


A flexibly mounted aircraft-propeller — power-plant combination can, 
under certain conditions, develop a dynamic instability in which the 
propeller axis together with the power plant precesses in a diverging 
whirl-type motion. The existence of such a phenomenon was recognized 
and briefly investigated by Taylor and Browne as early as 1938 (ref. l), 
but until recently has not been encountered in actual airplane configura- 
tions. An experimental investigation (unpublished) of this phenomenon 
was conducted on an aeroelastic model of an aircraft in the transonic 
dynamics tunnel at the Langley Research Center. These tests established 
the boundaries between stable and unstable operation of the system over a 
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range of such parameters as stiffness and damping of the engine mount, 
propeller speed, fuel loading, and so forth. 

The aim of this paper is to present a parallel analytical investiga- 
tion to shed further light on the basic mechanism of propeller whirl. 

The mathematical model treated represents a gross simplification of 
power-plant installations found in practice in that an isolated power- 
plant system is assumed to be flexibly mounted to a rigid backup struc- 
ture and the only aerodynamic forces considered are those acting on the 

propeller. It is hoped that, by following a tractable solution of the 
propeller-precession problem, reduced to its most elementary form, impor- 
tant nondimensional parameters will be identified and suitable design 
criteria will evolve by which the stability characteristics of arbitrary 
propeller — power-plant configurations can be readily assessed. 


SYMBOIS 


a i, b i 

b 

c e 

c l a 
C(k) 

Cm 


C n 

Cip 


C Y 


C Z 


aerodynamic coefficients in equation (19) 
propeller blade section chord 
viscous damping coefficient 

section lift-curve slope of propeller blade element 
Theodor sen's oscillatory lift function 

% 

pitching -moment coefficient, — 

pV^S'R 

Mz 

yawing -moment coefficient, 


pV 2 S'R 


T 

thrust coefficient, — — 


Fy 

side-force coefficient, 

|pV 2 S' 


vertical-force coefficient, 


ipV 2 S ' 
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!y^Z 


J 

k 

k e 
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M 

Me 


dissipation function (see eq. (5)) 

1 X Q 


angular momentum ratio, 


1 Y a B 


generalized aerodynamic force functions (see eq. (l8)) 

real and imaginary parts, respectively, of C(k) 
aerodynamic forces along Y- and Z-axes, respectively 

structural damping coefficient 

structural damping in pitch required for neutral stability 
of system 

S'ylf 

ratio of damping in yaw to damping in pitch, — 

60 


moment of inertia ratio, 




moment of inertia of propeller about axis of rotation 

moment of inertia of propeller — power-plant system about 
the pitch and yaw axes, respectively. 


propeller advance ratio. 


2hR 


reduced frequency of propeller blade element, ^ 

CJDaR 

reduced frequency of propeller — power-plant system, — 

distance from propeller plane of rotation to gimbal axes 

free-stream Mach number 

effective Mach number for propeller 
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My,M Z moments about Y- and Z-axes, respectively; positive sense 

is according to right-hand rule 

n propeller rotational speed, revolutions per second 

q^ generalized coordinate in Lagrange ! s equation 



generalized forces in Lagrange* s equation where q-j_ = 0 
or q ± = f 


R 

S* 

S 

t 

T 

U 

V 


propeller radius 
propeller disk area 

rotational stiffness of power-plant mount 
time 

kinetic energy of system or thrust 

potential energy of system or helical velocity of propeller 
blade element 

forward velocity 


X a 

Ws 

X,Y,Z 

0 


Po 

7 


equivalent air speed, > — — V 

\| P SL 

auxiliary coordinate in Euler axis transformation 
orthogonal space-fixed axis system (see fig. 2) 

orthogonal body-fixed axis system (see fig. 2) 

geometric blade angle measured from plane of propeller 
rotation 

blade angle to zero lift measured from plane of propeller 
rotation 

square root of ratio of yaw stiffness to pitch stiffness, 
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50, 6\|f 
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0,t 


K 


A 

P 

a 


T 


0) 

0)0 


"V 


0^,(1^, CDZ 


ft 


lag angle of oscillatory aerodynamic forces, tan - -*- 


-G(k) 


F(k) 


virtual work 

virtual displacement in 8 and a)/ 1 , respectively 
ratio of viscous damping to critical damping 
pitch and yaw angles, respectively (see fig. 2) 
effective pitch and yaw angles, respectively (see eq. (ll)) 

ratio of air to structural moment of inertia, — 


frequency ratio, — 
“0 


air density 

8 b 

propeller solidity at O. 75 R, 

nondimens ional distance traveled in propeller radii, — 

R 

frequency 


undamped natural frequency in pitch (at zero airspeed), 

\§i 

pY 

undamped natural frequency in yaw (at zero airspeed), 



components of total angular velocity of system about the 
X-, Y-, and Z-axes, respectively 

propeller rotational speed, radians/sec 
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Subscripts: 

av 

1,2 

o 

P 

q 

r 


average value 

forward and backward whirl modes, respectively 
amplitude of oscillation (eq. (20)) 
axes in plane of propeller rotation 

differentiation with respect to 

differentiation with respect to 


0V 

R 

jrV 

R 


rms root mean square 

SL sea level 

0.75R three-quarter propeller radius 

8,A|r pitch and yaw directions, respectively 


Dots over symbols indicate differentiation with respect to t; 
primes denote differentiation with respect to t. Partial derivatives 

dC Y dCy dc m 

are denoted thusly: Cy - - — , Cy = - — , C m = - — , and so forth. 

t 0 d0 1 dq 


DYNAMIC SYSTEM 


General Remarks 

The problem to be considered in this paper is the dynamic aeroelastic 
stability of a rigid engine -propeller combination which is flexibly mounted 
so as to permit angular deflections of the propeller axis in a vertical 
(pitch) plane and a horizontal (yaw) plane. (See fig. 1.) The primary 
forces and moments to be considered are those associated with small ampli- 
tude pitching and yawing motions of the system. These forces and moments 
are functions of the inertia of the system, the elastic and damping prop- 
erties of the mount, the gyroscopic moments due to propeller rotation, 
and the aerodynamic forces acting on the propeller. Thus, in contrast 
with actual airplane configurations, the wing structure on which the 


-si o ro h tr 1 
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power plant is mounted is assumed to be rigid and aerodynamic inter- 
ference effects induced by the wing and adjacent propellers are ignored. 

Consider first the free vibratory motions of the system when damping 
and aerodynamic forces are neglected. Assume that with a nonrotating 
propeller there is no coupling between the pitch and yaw degrees of free- 
dom; that is, natural vibrations in pure pitch and pure yaw can occur 
independently. When the propeller is spinning, however, gyroscopic moments 
are produced wherein an angular pitching velocity causes a moment about 
the yaw axis and vice versa. Thus, the significant characteristic of 
this system is that the pitch and yaw modes are coupled by the action of 
gyroscopic moments. (The coupling of natural modes of aircraft structures 
due to propeller gyroscopic forces is investigated by Scanlan and Truman 
in ref. 2 .) Gyroscopic coupling causes the natural frequency of the 
higher frequency uncoupled mode to increase and that of the lower fre- 
quency uncoupled mode to decrease. These modes are characterized by a 
precession-type motion. In the lower frequency mode the propeller axis 
precesses in a sense opposite to that of the propeller rotation and in 
the higher frequency mode the direction of precession and propeller 
rotation are the same. (See, for example, ref. 3*) In the present paper 
these modes will be referred to as the backward and forward whirl modes, 
respectively. 

Although gyroscopic effects tend to couple the pitch and yaw degrees 
of freedom, gyroscopic action alone cannot lead to an oscillatory 
divergent -type instability because the mechanism for adding energy to 
the system is lacking. The mechanism for an energy transfer, either into 
or out of the system, can, however, be found in the aerodynamic forces 
acting on the propeller. When the propeller axis is deflected in pitch, 
an aerodynamic vertical force and a yawing moment about an axis in the 
plane of the propeller are developed which, for small deflections, are 
proportional to the pitch angle. In addition to these static forces, 
other aerodynamic forces and moments proportional to the rate of change 
of the angular deflection are also present. Some of these air forces 
drive and others resist motion of the previously discussed whirl modes. 

A condition of neutral stability, that is, a self-sustaining oscillation, 
is achieved when the energy input due to aerodynamic driving forces on 
the propeller is exactly balanced by the energy dissipated by internal 
damping In the system. 

In the following analysis a modified Euler axis system together with 
Lagrange f s formulation of the differential equations of motion are used. 

The equations of motion are cast in a nondimens ional form from which 
approximate closed form expressions for conditions leading to neutral 
stability are derived and presented graphically. 
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Equations of Motion 

The coordinate system to he used in describing the motion of the 
propeller system is illustrated in figure 2. Two sets of orthogonal 
axes are considered. The X s -, Y s -, and Z s -axes are fixed in space with 

the origin at the gimbal center and the X— , Y— , and Z— axes pitch and yaw 
with the power plant but do not rotate with the propeller. The Euler 
angles relating the body system to space-fixed coordinates are the pitch 
angle 8 and the yaw angle \|r . The order of rotation in transforming 
from the space-fixed to the body-fixed coordinate system is as follows* 

L 

(a) Pitch through an angle 0 about the space-fixed Y s -axis to give 1 

X a , Y s , and Z coordinates. ^ 

9 

rj 

(b) Yaw through an angle i|r about the body-fixed Z-axis to give 
X , Y, and Z coordinates. 

It can be seen that the Z-axis is constrained to lie in the X s ,Z g plane y 

that is, the vertical plane. Positive sense is in accordance with the 
right-hand rule and is indicated by the direction of arrows in the fig- 
ure. The total angular velocity of the system is comprised of three 
vectors: a pitch rate 0, a yaw rate and the propeller rotational 

speed Q as indicated in figure 2. These vectors can be resolved into 
components a#, a*, and o£ along the X-, Y-, and Z-axes, respectively. 

The total kinetic of the system may be expressed 

T = | I x a X 2 + | I Y<^r 2 + | iz“z 2 (1) 

where I x is the polar moment of inertia of the propeller, and Iy and 
l£ are moments of inertia of the complete system about the Y- and Z-axes, 
respectively. In the present analysis it is assumed that I z = Iy. The 
kinetic energy of rotating engine parts could, have also been included in 
equation (l); however, since this energy is generally small, it has been 
neglected in the present analysis. The relations between the angular 
velocities in the body-fixed coordinate system and those in the space- 
fixed coordinate system are for small pitch and yaw angles (see fig. 2) 

ojy = fi + 8 sin ^ ** ft + 0^ 
coy = 8 cos \}r 0 

= ♦ 




-nI\d ro h 
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Q 


Thus, with equation (2), the expression for the total kinetic energy of 
the system for small pitch and yaw deflections about the gimbal axes 
becomes 


T 


= | l x (fi 2 + 2(10*) 



( 3 ) 


The potential energy associated with pitching and yawing deflec- 
tions is 


u = \ s 0 e 2 + i s t * 2 


(O 


where, as indicated in figure 1, Sq and are the rotational spring 

constants for pitch and yaw, respectively. 

Assume for the present that the damping in the engine mount system 
is defined by the structural damping concept commonly employed in 
flutter analysis. This type of damping produces a force in phase with 
the velocity of oscillation but proportional to the elastic restoring 
force. (See ref. 4 .) A dissipation function applicable to structural 
damping with sinusoidal motion may be written 


D = i s ese o2 i ^2 

2 cw 2 a> v 


( 5 ) 


where a) is the frequency of oscillation, and gQ and g^ are the 
damping coefficients for the uncoupled pitch and yaw modes, respectively. 

Although the choice of assuming a structural type of damping over, 
for example, viscous damping (force proportional to velocity), is at 
present arbitrary, since it is not known which is the better representa- 
tion of the damping in the system, it will be shown later that the type 
of damping assumed does have a significant effect on the stability 
boundary. 

The generalized forces are determined from the work done by the aero- 
dynamic forces acting through a virtual displacement of each generalized 
coordinate. In terms of the aerodynamic moments about the gimbal axes 
the total work done is 


6W = My&9 + 


( 6 ) 
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from which the generalized moments axe determined to be 


and 




% 

M Z 


( 7 ) 


where My and M Z are the propeller aerodynamic pitch and yaw moments 
about the gimbal axes. 


Now, the Lagrange equation is 


a for \ _ St_ du_ dD_ = 

dt fo q i/ Sq i Sq i ^ q i Qqi 

where is the generalized coordinate 0 or ^r. With the expressions 

for T, U, and D given by equations (3), (4), and (5), the differential 
equations of motion become 


Iy 0 + Iy N\(f + Sq 0 


i Y \|J - i x ne + s^\|r 


+ S 080 
0) 

+ CO 


0 =My I 

t = Mz 


(9) 


In order to determine My and M z in equation ( 9 ), it is convenient 

to consider next the aerodynamic forces and moments that act on the pro- 
peller. As stated previously, when the axis of propeller rotation is 
inclined relative to the direction of the free stream, aerodynamic forces 
and moments are generated which, for small deflections, are proportional 
to the angle of inclination and rate of change of the inclination angle. 

It is convenient to express these aerodynamic forces and moments in terms 
of nondimensional propeller derivatives which can either be determined 
experimentally or calculated by methods such as that given by Ribner in 
reference 5- On the basis of Ribner' s work the forces and the moments 
about axes in the propeller plane may be expressed as 
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f z = |^ 2 s’(c Z0 e + c^t + c Zr ^ 


My, P = pV 2 s-R(c m ^ + c mq 

r Y ’ + Cj/ + C Yq f) 


M Z , p - pV 2 S'R(c n0 e + C nr ^1 


( 10 ) 


J 


The derivatives C 7 and Cv„ are not considered in reference 5 tut are 
* J \|r H 

included here because these derivatives have been observed experimentally 
and can be shown to have a significant stabilizing effect on propeller 
whirl. The angles 9 and i|r are the effective pitch and yaw angles 
between the propeller axis and the relative wind. They consist of two 
parts, a geometric angle and an angle due to translational velocity of 
the propeller hub . Thus , 





J 


( 11 ) 


where Z is the distance from the rotation axes (gimbal axes) to the 
plane of the propeller. 

Because of symmetry the following reciprocal relations hold between 
the derivatives for pitching and yawing motions: 
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c z e = 

c ne - - c irty 

2y e - > d 2 ) 



With equation (lO) the moments referred to the gimbal axes became 

"N 

% = %,p ~ ZF Z 

and > ( 13 ) 

M Z - M z>p + ZF y I 

Before seeking solutions to these equations of motion it is advan- 
tageous to write equation (9) in terms of certain nondimens ional param- 
eters. For this purpose define a new independent variable t as the 
distance traveled in propeller radii 


T 


vt 

~R~ 


(i4) 


Then, in terms of t, the derivatives in equation (9) become 


) _ V d( 

R dT 


1 = X( )' 

IT 


( 15 ) 



i — I CVI CTs t — 


13 


and 


( 


) = 


(if 

W dT 2 



(16) 


With these substitutions and some algebraic manipulations the equations 
of motion may be written in the following nondimens ional form: 


e M + S ^ + g0 _i e' + k e 2 0 = nf e 

> 

V - | 9’ + ge ^ 7%' + 7 2 k 0 2 ^ = 


(17) 


The nondimens ional parameters in equation (17) are defined as 
follows : 

H ratio of polar moment of inertia of propeller to pitching moment 

of inertia of entire system, 

J propeller advance ratio, V/2nR 

reduced frequency, cdqR/V 

7 ^ ratio of yaw stiffness to pitch stiffn 
G ratio of yaw damping to pitch damping, 

k ratio of mass moment of inertia of cylinder of air to pitching 

moment of inertia of the structure, irpR^/jy 

A frequency ratio, o)/a>Q 




The generalized force functions are in the form 


fg = aQ0 + a 1 9' + a 20 " + bo^ + + l3 2'l f '' 


f. = -b 0 9 - b-^0' - b20" + + a i 1 t r * + a-gi" 


(18) 


where the a and b values in equation (l8) axe functions only of the 
propeller derivatives and the gimbal axis location. These coefficients 
are: 


lip 

0 ~ " 2 H CZe 
a l - Cm q + c Ze 


a 2 


b 0 


- 1 % 


( 19 ) 


= - c z 

R\2 


+ b 0) 


b 2 





It should be pointed out that the present analysis does not consider 
the apparent mass effects of air surrounding the propeller, which, if 
accounted for, would lead to additional terms in the a 2 and b 2 

coefficients . 



— ^ vo ro m tr* 




Evaluation of Propeller Aerodynamic Derivatives 

The aerodynamic derivatives used in this analysis were calculated 
for a propeller having the characteristics shown in figure 3. All 
derivatives were evaluated by Ribner’s method (ref. 5) with the excep- 
tion of which is not treated by Ribner. Such a derivative has, 

however, been observed experimentally and is believed, because of its 
stabilizing effect, to be of sufficient importance in the analysis of 
propeller whirl to merit some consideration. An approximate derivation 
of the derivative is presented in the appendix. 

Because, as will be shown later, the thrust developed by a pro- 
peller has very minor effects on the aerodynamic derivatives of interest, 
wind-milling propellers were used here as they were in the experimental 
wind-tunnel tests (unpublished) . From experimental data, the relation 
between the propeller advance ratio J and the geometric blade angle 
P 0 was determined and is plotted in figure 4. Figures 3 and 4 

together with data in table I provide the information needed to cal- 
culate propeller derivatives for the windmilling case. 

The calculated propeller derivatives plotted against |3q y^p are 

presented in figure 5* Also shown for comparison are the derivatives 
calculated for a thrusting propeller for Pq y^p = 46°. 

In order to account for the effect of Mach number on the coeffi- 
cients in figure 5* an approximate Mach number correction factor 

1 

- — ■ — , calculated by the method of reference 5 ; is shown plotted 

\l 1 - M e 2 

in figure 6 against P 0 . 75 R f*°r various free-stream Mach numbers. 

(The symbol M e represents an effective Mach number on the propeller.) 

This correction is to be applied as a multiplying factor to the deriv- 
atives determined for the incompressible case. 
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DYNAMIC STABILITY CONSIDERATIONS 


Conditions for Neutral Stability 

In order to find the conditions for which the system becomes neu- 
trally stable, sinusoidal motion is specified for 8 and \Jr in the 
equations of motion. Thus, when 


0(t) = 0 o e lajt 


or in nondimensional form 


0 ( t ) = e o e iAk 0 T 

and 




J 


( 20 ) 


are substituted into equation (17), & nontrivial solution is obtained 
when the determinant of the coefficients 


[y(i - * 2 ) + K(* 2 k e 2 a 2 - ao)J + ifg 0 k 8 2 - tXkga^j - b 0 )J + ijxkg^ - *b^ 

^b 0 - k 2 k e 2 b 2 j] + i Ak^Kb-L - 2^ [^e 2 ( y 2 - * 2 ) + K(x 2 k 9 2 a 2 - Sq)] + i [ggGkg^ 2 . nAkga^J 


( 21 ) 


is set equal to zero* 

When equation (21) is expanded and the resulting real and imaginary 
parts equated to zero, it is found that in the real equation gg always 

appears either multiplied by itself or by one of the aerodynamic terms* 
By neglecting these second-order terms, which appears to be a valid 

assumption for the range of k- and g-values of interest, the real part 

o 

of the expanded determinant becomes a quadratic in A . The two roots 

p 

of ~K , corresponding to the forward and backward whirl modes, are, 
respectively. 
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Ai 2 = 


-B + l/i 2 - 4AC 
2A 


A 2 2 = 


■B - j/j 2 - 4AC 
2A 


y 


( 22 ) 


where 


A 



2 k 3*2 t 



+ hr 


) 


B = kr 


-(1 + 7 2 Ae 2 - (f) 


+ 2Ka, 


■0 


+ k c 


2bp j + a 2^ + 


- 2k 


( a 0 a 2 + b 0 b 2 


K (a]_2 + b]_2^ 


C = j 2 K 


Ka r 


1 + y 2 )\ 


K ( a o 2 + b o 2 ) 


The terms beneath the bars are, for conditions of the present analysis, 
considered to be negligible and have been ignored. 

The amount of damping required for neutral stability, is found by 
equating the imaginary part of the expanded determinant to zero and 
solving for gQ 


g 


e 


Xk 


a l 

to 2 2 2 2^ 

V - A + 7 k 9 ] 

+ 2 i 

(b 0 - vV’ 

- 2lt (^cAi + b 0 b : 

) + 2K^a^a2 + 


k e 2 (r 2 - > 2 ) + k e 

2 (l - A 2 c 

+ K 

(*A 2 - a 0 ) 

\(l + 7 2 o) 


(23) 


The ~ superscript is used to denote a neutral stability condition. 

The procedure for determining stability boundaries from equa- 
tions (22) and (25) is as follows: For a given set of parameters 

equation (22) is first evaluated to give the forward and backward whirl 
frequencies ?\p and Only positive roots are used in order to be 
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consistent with the motion assumed in equation (20). With these 
frequencies substituted separately into equation ( 23 ) the damping 
required for neutral stability is determined for each whirl mode. 

Figure 7 shows typical plots of A-, 0 and g n against the 

0^2 

parameter yRjQQ for a system having symmetrical damping (G = 1. 0) 

and asymmetrical stiffness (7 = 1.4). Other parameters used are pre- 
sented in figure 5 and table I. If for a given set of conditions the 
system damping is less than that indicated by the appropriate g^ 

stability boundary in figure 7, the system will be unstable. Since 
the damping in an actual physical system is always positive, the fig- 
ure shows that propeller whirl, if encountered, will stem from the 
backward whirl mode. 


In order to gain a better understanding of the significance of 

various terras in the stability boundary equation, it is instructive 

to make further approximations to arrive at a considerably simpler 

expression for g . For this purpose consider the special case of 
9 

equal stiffness and damping in the pitch and yaw directions, that is, 
7=1-0 and G = 1.0. If in the frequency equation (eq. ( 22 )) all 
aerodynamic terms are neglected and an angular momentum ratio, defined 
ivfi 

as E = - is assumed to be less than unity, the frequency relation 

1^9 

reduces to 


A 


1 


1 




« 1 



(24) 


When these relations for A are substituted into equation ( 23 ) and 
the aerodynamic terms of lesser importance dropped, the following 
approximate expressions for g Q are obtained 

Forward whirl: 



Backward whirl: 
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Where, as in equation (19), 



Figure 8 indicates that these approximate expressions for ^ 

and g agree very closely with the more exact equations (eqs. (22) 

y l,2 

and (23)) for the case under consideration, that is, a symmetrical 
distribution of stiffness and damping in the system. 

Note that for the aerodynamic coefficients and gimbal axes loca- 
tion used (see fig. 5 and table I) the term a^ I s always negative 

and the term bQ is always positive. Therefore, equations (25) and (26) 
show that for these conditions the damping required for neutral sta- 
bility can become positive only for the backward mode. This result 
indicates that, if an instability occurs, it must develop in the back- 
ward whirl mode, a conclusion consistent with wind-tunnel results and 
the theory in reference 1 . 


Subcritical and Supercritical Response 

The equation for gg defines only the boundary between stable 

and unstable operation of the system and does not necessarily indicate 

the degree of stability or the severity of an instability when the 

actual damping is other than g Q . 

y 


As a qualitative indication of the behavior of the system for 
conditions on either side of the boundary, figure 9 illustrates some 
typical transient motions of the system for three values of damping. 
These curves, obtained on an analog computer, represent the paths 
traced by the propeller hub after being released from an initial lateral 
displacement. The cases shown are for equal stiffness and damping in 
the pitch and yaw directions. Note that the oscillatory motion is 


damped when g Q = 0. 08, neutrally stable when 
y 

and divergent when gQ = 0 . Ok. 


ee = Sq 2 = °- o6 > 


The path described by the propeller hub may be considered as the 
superposition of two rotating vectors. One vector rotates about the 
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origin of the Y S ,Z S coordinates in a clockwise direction at the fre- 
quency of the forward whirl mode and the other vector rotates in a 
counterclockwise direction at a lower frequency corresponding to the 
backward whirl mode. For the particular initial condition used in fig- 
ure 9 a corner occurs in the path each time the two rotating vectors 
point in the same direction. Note that in each case the backward whirl 
(counterclockwise) is the predominant mode and that the forward mode is 
heavily damped. 


Generalized Stability Boundary Plots 

Description. - Figure 10 gives the effects of varying some of the 
nondimens ional parameters considered in this analysis. These charts 
show gQ^ plotted against the parameter V/Rug for a range of values 

J > 7 > and G. Only gg values for the backward whirl mode are 

presented because in all cases considered this was the mode in which 
the instability occurred. In these calculations the following condi- 
tions were assumed: windmilling propeller, incompressible flow, fixed 

pitch, and yaw rotation axes at the position of the model gimbal axes 
(Z/R = 0 . 3778 ) , and inertia parameters corresponding to the system at 
a flight altitude of 15,000 feet standard atmosphere. (See fig. 5 and 
table I for numerical values. ) 


Before utilizing these nondimensional plots it is of interest to 
consider trends that might be expected when conditions are different 
from those assumed in figure 10. Specifically, the effects of varying 
the air density, Mach number, propeller thrust, gimbal axes location, 
and the assumed nature of the damping in the engine mounts are discussed. 


_ Air density .- First consider the effect of air density on the sta- 
bility boundary. For the case of symmetrical stiffness and damping it 
is seen from the approximate relation given by equation (26) t ha t g 0 

is simply proportional to the nondimensional moment of inertia param- 


eter 


k = 


rtpR' 


Thus, in order to determine approximately the sta- 


bility boundary for altitudes other than that assumed in the figure, 
the value of read from figure 10 can be multiplied by the ratio 


of the air density at the desired altitude to the density used in this 
analysis (15,000 feet). 


Mach number .- The effect of Mach number on the stability boundary 
is indicated in figure 11 for J = 2 . 6 , G = 1 . 0 , and y = 1 . 0 . These 
results were obtained by applying the Mach number correction shown in 
figure 6 to the propeller aerodynamic derivatives. Note that for a 
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given mount damping and stiffness (stiffness is proportional to oq 0 2 ) 

the effect of increasing Mach number is to reduce the velocity for which 
the system becomes unstable. The higher the Mach number the more sensi- 
tive the system becomes to this unfavorable effect. 

Propeller thrust .- Shown in figure 12 is a comparison of the cal- 
culated stability boundary for a thrusting and a windmilling propeller. 

The propeller coefficients used for the thrusting case are those plotted 
in figure 5 and correspond to cruising thrust at 15,000 feet altitude 
with an advance ratio J = 2 . 22 . Note that the differences between these 
two stability boundaries are hardly discernible in the figure. Conse- 
quently, there appears to be good justification for using windmilling 
in place of thrusting propellers in the wind-tunnel tests (unpublished). 

Gimbal axes location .- The stability boundaries shown in figure 10 
are for gimbal axes located 0. 3778 propeller radius behind the propeller 
plane of rotation. In order to examine the sensitivity of these boundaries 
to gimbal axes location, the parameter l / R was varied in figure 13 over 
a range from 0 to 0 . 80 . In these calculations it was assumed that only 
the aerodynamics were affected by changes in z/R, the moments of inertia 
remaining fixed. These results indicate that increasing the distance 
between the propeller plane and the gimbal axes has a stabilizing effect. 
Note also from equations (25) and (26) that as this distance is increased 
beyond about 2 propeller radii the sign of bQ changes from positive to 

negative. In this case, whirl instability would be expected to develop 
in the forward whirl mode; however, limited calculations indicate that 
before reaching this condition, which occurs at very large values of 
V/Rcb e , the system would probably suffer a static aeroelastic divergence. 


Damping concept . - As pointed out earlier, the type of damping con- 
sidered for the system produces a force proportional to the amplitude 
of displacement but independent of the velocity of motion. This concept 
of damping has been found to give a better description of damping meas- 
ured on aircraft structures than does viscous damping which produces a 
damping force proportional to velocity. It is of interest, however, to 
compare the stability boundaries determined for these two familiar types 
of damping. For a sinusoidal pitching motion of frequency o>, It can be 
shown that, for a given amplitude of oscillation, structural damping and 
viscous damping dissipate the same energy when 


s 0 


2t a_ 


(27) 


where 
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i 


0 


2I tf°0 


is the viscous damping coefficient relative to critical damping. In 
this equation cdq is the undamped natural pitching frequency of the 
system with a nonrotating propeller and co represents the damped natural 
frequency. If the propeller is not rotating and the system is lightly 
damped, o> =» cjdq and therefore from equation (27) structural and viscous 

damping are simply related as gQ ~ 2£g. When the propeller is spinning, 

however, because of gyroscopic coupling, co may be considerably dif- 
ferent from <x> Q , in which case the damping forces will depend on whether 

the damping is of a viscous or of a structural nature. 

An indication is given in figure Ik of the shift in the stability 
boundary due to using a viscous rather than a structural description 
of the damping in the engine mount. The viscous damping curve was com- 
puted by replacing g Q in equation (9) by the equivalent expression 

given in equation (27) (similarly for g^). As is evident in the fig- 
ure, viscous damping is not as effective as structural damping in the 
prevention of instability since the unstable mode, being the backward 
whirl mode, has a frequency less than that for the case of a nonrotating 
propeller. Therefore, to absorb a given amount of energy, the viscous 
damping coefficient 2 £q must be greater than the structural damping 

coefficient gQ. 


Effects of Unsymmetrical Stiffness and Damping 


Experimental results for the stability boundaries of systems having 
different relative lateral and vertical stiffnesses are sometimes cor- 
related by an effective stiffness defined as the root mean square of the 
two stiffnesses 


S 


rms 



(28) 


Thus, within the accuracy of this criterion, the stability char- 
acteristics of systems having arbitrary relative stiffness in pitch 
and yaw can be inferred from the behavior of a symmetrical system having 
the same root -me an- square stiffness. It is of interest to determine the 
scatter that might be expected when using this method of correlation. 
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These results would provide an indication of the range of stiffness 
ratios over which the previously developed simplified stability boundary 
equation for a symmetrical system (eq. ( 26 )) would be applicable. 

Figure 15 (a) presents a family of stability boundaries showing the 
variation of root -me an- square stiffness with equivalent airspeed for 
various ratios of yawing stiffness to pitching stiffness. A structural 
damping of 0.03 was assumed in both pitch and yaw and a propeller speed 
of 1,020 revolutions per minute was used. These plots were determined 
graphically from the nondimensional working charts in figure 10. Fig- 
ure 15 (a) indicates that at the lower root-mean- square stiffness values, 
the root-mean-square stiffness criterion is a reasonable criterion to 
use in defining the stability boundary over the range of stiffness 
ratios shown. Also, the figure points up the interesting fact that, 
for a given root -mean- square stiffness, it is the symmetrical system 

= l.o] that has the lowest critical airspeed. Thus, for a given 
s e / 

root-mean- square stiffness, engine-mount stiffness asymmetries appear 
to produce a beneficial effect; however, the range of parameters for 
which this conclusion applies has not been established. 

The structural damping in figure 15(a) was assumed to be equal for 
both modes. In order to illustrate the effect of unsymmetrical damping 
on the stability boundary, the average damping 


£>av 


ge + gfr 
2 


( 29 ) 


was held constant at 0.03 and the damping ratio G was varied from 0.5 
to 2.0. The results are presented in figure 15(b) which shows stability 
boundary plots of stiffness ratio against equivalent airspeed for a 
root-mean- square stiffness of 7.36 x 10 ^ and 12.3 x 10 ^ in-lb /radian. 

The G = 1.0 curve represents a cross plot of figure 15 (a). If the 
stiffness is symmetrical, the average damping is a gpod measure of the 
effective damping of the system since the stability boundary for each 

c» 

damping ratio shown intersects at -Jfc = 1.0. However, when the stiff - 

s e 

ness ratio is different from unity and the average damping held con- 
stant, it is seen that the critical speed is increased when the damping 
in the less stiff direction exceeds that in the more stiff direction. 
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COMPARISON OF THEORY WITH EXPERIMENT 


General Trends 

In the present sections, a comparison is made between the results 
of this analysis and experimentally determined propeller whirl char- 
acteristics. Although the analysis shows the same trends as the experi- 
mental data, close agreement is not to be expected since the theory 
considers only an isolated nacelle, and thus neglects wing motion and 
aerodynamic effects of the wing and adjacent propellers. It should also 
be noted that the actual structural damping present during the experi- 
mental tests could not be precisely determined because, in order to 
obtain consistent results, it was necessary to measure damping with the 
nacelle system mounted to a rigid backstop rather than to the flexible 
wing structure. Also aerodynamic drag loads on the gimbal bearing may 
tend to increase the damping above that measured in still air. 


Whirl Frequency 

A comparison of calculated and measured whirl frequencies against 
propeller speed is shown in figure 16 for two stiffness ratios. Aero- 
dynamic forces have a small effect on the frequency of whirl and were 
neglected in determining the theoretical curves shown. The figure 
clearly indicates that the measured whirl frequencies stem from the 
backward whirl mode. 


Propeller-Speed — Airspeed Stability Boundary 

Figure 17 shows the measured and calculated effects of propeller 

speed on the stability boundary. The stiffness is 8.09 X 10^ inch-pounds 
per radian and the structural damping is approximately g Q 0.014. The 

two experimental stability boundaries are for the number k (outboard) 
nacelle and the number 3 (inboard) nacelle. Thus, the three curves in 
the figure represent three different levels of flexibility in the nacelle 
mount backup structure. Nacelle 4, being located outboard on the wing, 
has a less rigid backup structure than nacelle 3 . The theoretical curve 
is, of course, a limiting case wherein the backup structure is completely 
rigid. From these curves, it may be concluded that the effect of 
increasing wing flexibility in this case is to increase the airspeed for 
which the system becomes neutrally stable. Also, both theory and exper- 
iment show that an increase in propeller speed is accompanied by a reduc- 
tion in critical airspeed. 


B 
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Damping-Airspeed Stability Boundary 


The variation of structural damping required for neutral stability 
with airspeed is shown in figure 18 for conditions of a propeller speed 

of 1,020 revolutions per minute and a stiffness of 8.09 X 10^ inch-pounds 
per radian. If, for a given airspeed, the structural damping is greater 
than that indicated by the stability boundary, disturbed motions of this 
system will dampen out. If, on the other hand, the damping of the mount 
falls below the stability boundary curve, disturbed motions will develop 
into an oscillatory divergence in the backward whirl mode. Here, as in 
figure 17 , the same indications regarding the effects of wing flexibility 
are evident. Note that, for low values of structural damping, the slope 
of the stability boundary curve is small (see fig. 10 ) and therefore, 
in this region, the critical airspeed is very sensitive to changes in 
damping. 


Stiffness-Damping Stability Boundary 

As a final comparison of theory (rigid backup) and experiment 
(unpublished), figure 19 shows the variation of root-mean- square stiffness 
with damping required for neutral stability at V e = 504 knots. The 

theoretical curve is for a propeller speed of 1,020 revolutions per 
minute and symmetrical stiffness and damping. The experimental results, 
are for the number 4 nacelle and fall within the band indicated by the 
crosshatched area. The apparent scatter in the experimental results 
is to be expected since the data shown represent a composite of all 
data points irrespective of propeller speed or relative stiffness and 
damping in the pitch and yaw directions. The theoretical results are 
again conservative in that for a given damping a higher stiffness is 
required for stability by the theory than was determined experimentally. 


CONCLUDING REMARKS 


An analytical investigation has been made of the dynamic stability 
of an aircraft engine-propeller combination flexibly mounted to a rigid 
backup structure. The conditions for which the system undergoes self- 
excited precession (propeller whirl) are examined. Results are pre- 
sented in the form of nondimensional plots from which stability bound- 
aries for various combinations of stiffness of the power-plant mount, 
structural damping, and propeller speed can be readily determined. Sta- 
bility is found to depend strongly on the damping and stiffness in the 
system. It is found that the theoretical results, which do not account 
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■^or the effects of wing response, show the same trends as observed in 
wind-tunnel tests of an aeroelastic model but indicate instability for 
lower velocities than do the model tests. Thus, wing flexibility appears 
to increase the damping of the system. ~ ' ' 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Field, Va., November 14 , i960. 
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APPENDIX 

PROPELLER VERTICAL FORCE DUE TO YAW DERIVATIVE C z ^ 


The existence of the derivative C 7 can be attributed to phase 

lags of aerodynamic forces acting on the propeller. When the propeller 
axis is deflected relative to the free- stream direction, as in yaw, the 
angle of attack of a blade element varies harmonically at a fundamental 
frequency corresponding to that of the propeller rotational speed. 
Because of an aerodynamic phase lag associated with these angle-of- 
attack oscillations, the total force vector in the plane of rotation 
due to a yaw angle has a small vertical component as well as the usual 
horizontal component. Thus, if the side force due to yaw undergoes a 
lag 5, the vertical force due to yaw may be expressed as 

CZ-,!; ~ Cy tan 5 (Al) 


A rough estimate of the magnitude of 5 can be obtained by simply 
modifying the slope of the section lift curve of the propeller blade 
element by Theodorsen f s C(k) function (ref. 4) for oscillatory flow 

C(k) - F(k) + iG(k) (A2) 

The lag angle 6 is given by the expression 


6 


tan“l 


-G(k)~ 
F(k) _ 


(A3) 


This expression for 5 must, at best, be viewed as a very approximate 
one since C(k) is derived on the assumption that the shed vortices be 
in a plane whereas in the case of a propeller the wake is shed in a 
helical pattern. 


The real and imaginary parts of C(k) are functions of a reduced 

frequency k = which, of course, varies along the propeller radius. 

For the present purpose assume an effective reduced frequency for the 
propeller to be that of the blade element at 0. 75R« It can be shown 
that the velocity of this blade element along a helical path expressed 
in terms of the rotational speed and the advance ratio is 


U 0. 75R 



(AM 
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Hence, the effective reduced frequency becomes 



From tabulated values of C(k) the phase lag and thus C 7 


(eq. (Al)) 


can be determined, for specified values of J. It is found that the phase 
lag varies from 11.9° to 9- 4° over a J range of from 1.8 to 4-. 2. 
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TABLE I 

SUMMARY OF NACELLE- PROPELLER CHARACTERISTICS USED IN ANALYSIS 


1^, slug- ft 2 
Iy> slug- ft 2 
, ft ... 

0.75R' ft * 

P , slug/ f 

z/r .... 

a 



Spinner- nacelle fineness ratio 
Spinner radius/Propeller radius 


175 
l, 575 
6.75 
1.458 

0. 001496 
0. 3778* 
0. 1834 

2rt 

6.0 

0. 181** 


*In figure 13 Z/R is varied. 

**Values assumed for use in reference 5. 




Propeller twist distribution, deg 


•6 


.h\— 




Figure 3*“ Distribution of propeller twist angle and blade width for the 
airplane propeller of this analysis. Propeller dimensions: blade 

tip width, 1.458 ft: propeller radius, 6.75 ft: number of blades, 4. 
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of advance ratio with geometric blade angle at 
blade radius for windmilling propeller. 




Aerodynamic derivative 
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Figure 5*“ Variation of propeller aerodynamic derivatives with geometric 
blade angle at three-quarter blade radius, M = 0. 
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Figure 6.- Variation of Mach number correction factor for propeller 
aerodynamic derivatives with geometric blade angle at three-quarter 
blade radius. 
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Figure 7«- Typical plots of frequency and damping required for neutral 
stability for various propeller advance ratios. S^/Sq = 1 * 96 ; 

G = 1.0; 7 = 1.^. 
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Figure 10.- Concluded 












o l 2 3 h S 


v 

Rw e 

Figure 12.- A comparison of the stability boundaries for thrusting and 
windmill ing propellers. J = 2.22; G = 1.0; y = 1.0. 




Figure l4 . - Comparison of stability boundaries for viscous and structural 

damping. J = 2.6; G = 1 . 0 ; y = 1 . 0 . 
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Figure 15*- Concluded. 





Propeller speed, rpm 
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Figure 17 .- Stability boundary plots of propeller speed plotted against 
airspeed for propeller whirl. G = 1.0; 7 = 1 . 0 ; gQ^ 0.014; 

Sq = 8.09 X 10^ in-lb/radian. 




Figure 18 .- Stability boundary 


against airspeed for propel] 
S 0 = 8.09 X 10 6 in-lb/radiai 













